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From the many-body T-matrix the condition for a medium-dependent bound state and its binding 
energy is derived for a homogeneous interacting Bose gas. This condition provides the critical hne 
in the phase diagram in terms of the medium-dependent scattering length. Separating the Bose 
pole from the distribution function the influence of a Bose condensate is discussed and a thermal 
minimum of the critical scattering length is found. 
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The discovery of Feshbach resonances in gases of ultra- 
cold bosons in 1998 [l|, y, IsJ has provided an important 
tool to analyze interacting Bose systems. Near these reso- 
nances it is possible to tune the interaction and especially 
the free scattering length aq vifith an external magnetic 
field B. In the vicinity of the resonance dX B — Bq the 
scattering length is [J| 
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as shown in Fig. [T] Onr is the scattering length far away 
from the resonance and Ai? oc l/cnr describes the width 
of the resonance. A Fermi gas near a Feshbach resonance 
can be driven through a transition from a Bose-Einstein 
condensate (BEC) of two-particle bound states for ao > 
to a BCS state of Cooper pairs for ao < 5] . For bosons 
the influence of the interaction on the Bose condensation 
is of main interest. Interacting Bose gases at ultra- low 
temperatures are expected to consist of unbound, bound 
and condensed bosons, furthermore one expects an in- 
fluence of the interaction on the critical temperature and 
density of Bose condensation, see citations in [6[ . Here we 
focus on the formation of bound states in the presence or 
absence of a Bose condensate. We will derive the condi- 
tion for bound states in terms of the medium-dependent 
scattering length to discuss the phase diagram and the 
binding energy. 

The two-particle scattering is described with the many- 
body T-matrix in ladder and quasi-particle approxima- 
tion 
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The influence of the surrounding particles is represented 
by the distribution function /p. The total momentum Q 
reflects the center-of-mass motion of the scattering parti- 
cles relative to the medium, while q is their relative mo- 
mentum. The medium is assumed to be a homogeneous 




FIG. 1: Scattering length in the vicinity of a Feshbach reso- 
nance 



ideal Bose gas with a distribution [7| 
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where tiq is the condensate density and F is the total spin. 
Since we assume that the interaction does not change the 
spin there are no degradation factors in the T-matrix ([2]). 
In the normal state no = and Eq. (j4|) yields the chemical 
potential /x as a function of the temperature and density. 
In the superfluid state /i = and Eq. ^ determines 
the condensate density uq. At the critical point ^ = 
and no = and from ([4]) follow the critical temperature 
Tq and the critical density uq. The dependence of the 
critical properties on the interaction will be neglected. 

At low temperatures only s-wave scattering at small 
momenta is important, furthermore we want to concen- 
trate on bound states near the continuum edge. There- 
fore and for the sake of simplicity we neglect the range 
of the interaction and assume the interaction to be a 



contact interaction for which the potential is indepen- 
dent of the relative momenta, Vpq = V. Accordingly, the 
T-matrix is also independent of the relative momenta, 
Tpq{Q,oj) =T{Q,uj), and Eq. ^ simplifies to the al- 
gebraic relation T = V/(l — QV) with the two-particle 
progator 
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Obviously we can split Q into the free propagator Qq re- 
sulting in the / — > limit and the medium correction 
Qra oc /. The free propagator diverges and a cutoff is 
necessary. To circumvent this cutoff, we introduce the 
vacuum T-matrix To = V/(l — QoV), which determines 
the free scattering length 
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The strength of the interaction is now described by the 
free scattering length, which is interpreted as the relevant 
physical quantity tunable near the Feshbach resonance. 
Solving ^ yields 
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The inverse of the cutoff go is proportional to the range of 
the interaction, i.e., contact interaction means qo — > oo. 
Postulating a finite an, it is necessary to renormalize the 
interaction strength [g, |3, Isl, llO[ , such that 
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FIG. 2: Binding energy for the bound state of a liomogeneous 
Bose gas with contact interaction. Fixed parameters are given 
above the plots, nc = 0.06s"^, To = O.QOss, ac = 0.06s 
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follows from ([7]). Frequently the pseudopotential [4l.lll| 
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is used to describe the dependence of the interaction 
strength on the scattering length for contact interaction 
[12I ll3J . The difference to ^ is that the leading term 
of ^ with respect to a large cutoff qo is always nega- 
tive and independent of ao, i.e., the contact interaction 
is always attractive, while the pseudopotential is attrac- 
tive for ao < and repulsive for ao > 0. The reason for 
this difference is that the pseudopotential is only valid 
up to first-order Born approximation 
approximation one neglects Q, i.e., T 
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. Within this 
V, such that 

the diverging terms in the denominators of ([7]) and ([5]) 
vanish. 

The in-medium T-matrix can be expressed by the free 
one T = 7o/(l — ^m^)- Accordingly, the in- medium 



scattering length is 
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with ac > 0. We use here the definition of the many- 
body scattering length of lj| instead of the definition 
used in [y, Isj. Instead of a divergence the many-body 
scattering length for Q = therefore has a zero at the 
critical point of Bose condensation, i.e., ac — 0. Splitting 
([5|) in another way 
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FIG. 3: Phase diagram for the appearance of the bound state 
(BS) of a homogeneous Bose gas with contact interaction. 
The vertical lines mark the onset of Bose condensation for 
the corresponding parameters 



one obtains for the in-nicdium T-matrix 
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A bound state, i.e., a pole of the T-matrix, is therefore 
possible if 
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where the corresponding binding energy lob < is mea- 
sured relative to the continuum edge, i.e., J^ . The 
binding energy is shown in Fig. [21 The fixed parameters 
in the plots define the corresponding length scale s, i.e., 
Q ^ p^ ^ l/s, T ^ Es = h^/2ms'^ or n = 0.05s-^ re- 



spectively. According to pO]) the condition for the bound 
state a > is satisfied in two situations. In the first 
case ao > the interaction is strong enough to form the 
bound state already in the vacuum. In the second case 
ao < — ac the bound state is induced by the medium. 
For — ac < a < the interaction is attractive but insuffi- 
cient to form a bound state. For bosons near a Feshbach 
resonance this means that in addition to the bound state 
in the ao > region, a bound state is also possible on 
the other side of the resonance for ao < — ac- The ap- 
pearance of the medium-induced bound state is also sig- 
naled by the divergence of the medium-dependent scat- 
tering length pO|) at ao ~ — ac [15|. From these three 
cases the dependence of the bound-state region on the 
density, temperature, total momentum, scattering length 
and spin follows as shown in Fig. [31 
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FIG. 4: Inverse critical scattering length versus temperature 
for different approximations (upper plot), and for different 
total momenta for a Bose gas with condensate (lower plot) 



For large |ao| the binding energy converges towards a 
finite value wc as seen in Fig. [He), which follows from 
J{Q,uJc + h^Q^ l^m) = 1/ac- This convergence can be 
explained by the convergence of the medium-dependent 
scattering length (|10p and also by the convergence of the 
interaction strength ([5]) for large |ao|. The interaction 
strength ^ also shows that the interaction is stronger 
for positive ao and therefore JwbI is larger in this case. 
The region where a bound state is possible spreads with 
increasing density, as can be seen in Fig. [31a). Fig. [21a) 
shows that for increasing density \uj-q \ increases, too. This 
behavior can be explained by the increasing influence of 
many-body effects with increasing density. Due to the 
Bose enhancement the formation of bound states is sup- 
ported. On the other hand thermal fluctuations hinder 
the formation of bound states, which is shown by the 
shrinkage of the bound-state region in Fig. [31b) and the 
decrease of \oj-b\ va Fig. [2lb) with increasing tempera- 
ture. The motion of the scattered particles relative to the 
medium has a similar effect as Figs.[2lc) and[3lc) show. In 
the limit of vanishing density the bound-state condition 
is ao > and the binding energy is wbo = — h^ /ma^. 
This simple result can only be explained by the potential 



TABLE I: Binding energies for different species. The singlet (s) and triplet (t) scattering lengths ao are from [j if not marked 
differently. The length scale s is chosen so that T = £s = 0.5/iK fcB~43peV. The binding energy ub as follows from (|14p for 
Q = and n = 0.85nc is compared to its vacuum value ujbo = — fi? /ma^, ac ~ 64x 10~^s. The total spin F is that of the 
hyperfine state with lowest energy. 
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-(112.2 ±2.2) X 10"^ 
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-15.4 ± 1.2 
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no bound state 
-(19.6 ±1.7) X 10"^ 

-0.039 ± 0.022 
-1.44 ±0.11 
-15.3 ±1.2 
-21.2 ±0.6 
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30 ±6 [161 
23.8 ± 0.4 
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3.1 X 10"^ 
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-1800 ±500 
-2000 ± 400 
-1.04 X 10'^ 
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(j8]). If one would follow the philosophy of pseudopoten- 
tials instead, one would have repulsion for ao > 0, and no 
bound states. In Figs.[2]d) and[3]d) one sees that with in- 
creasing total spin JwbJ the bound-state region decreases. 
The reason is that with increasing total spin the density 
of states also increases and therefore the occupation of 
states for a given density decreases. The effect is there- 
fore similar to that of a decrease of density. 

A remarkable feature in Figs.[3]c) and d) is that a maxi- 
mum of the critical line appears at some temperature T^,^. 
This means that in these cases the region of bound states 
becomes smaller if the temperature is decreased further, 
which is in contrast to the behavior above Tex. To un- 
derstand this effect we observe from ^ and pil| that 
in the superfluid state 1/ac can be split into two parts, 
1/ac = l/ag'"<^ + l/ag'". The first part 1/ag'"'^ oc no/Q^ 
bears the contribution from the condensate. The second 
part l/a|f'' oc / represents the uncondensed Bose gas 
vanishing as the temperature approaches zero. While 
aff' diverges at T = and ag'"^ diverges at T = Tq the 
scattering length ac remains finite as Figs. [3]c) and d) 
show. If the momentum Q is large enough the extremum 
appears already above the critical temperature as seen in 
Fig. m The condensate part vanishes at Tq and the gas 
part is the only contribution above. For large total mo- 
mentum the Bose distribution is well approximated by 
the Boltzmann distribution. This allows us to calculate 
Qra explicitly which yields that the extremum appears at 
Tex ~ 0.22 h^Q^ /Am. In other words we have a mere 
thermal effect. When the mean motion characterized by 
the total momentum comes in resonance with the thermal 
motion we observe an extremum in the critical scattering 
length. 

The ladder and quasi-particle approximation proposed 
in this paper cannot describe the dynamic formation or 



breaking of bound states. This is due to the fact that 
for the formation or destruction of a bound state an ex- 
change of energy and momentum either with a third par- 
ticle or with the medium has to be allowed. Therefore 
either more diagrams or a self-consistent spectral func- 
tion would have to be included in the T-matrix. This is 
the reason why the pole of (|13|) lies on the real axis, i.e., 
has no imaginary part, and therefore the bound state has 
an infinite life time. On the other hand Figs.[2]a), b) and 
e) show that in the normal state |ljb | is of the order of the 
thermal energy, for ao < — ac and for oq S> ac, i.e., near 
the resonance. This makes the bound states very instable 
towards collisions with medium particles and therefore 
limits their life time. However experimental experience 
shows, that bound-state and cluster formation has even 
to be suppressed by decreasing the density to be able 
to directly investigate Bose condensation [17]. If Bose 
condensation shall be reached by decreasing the temper- 
ature, then bound states are always possible before Tc is 
reached for any interaction, as Fig. [3]c) shows. The only 
exception is the ideal gas case, i.e., ap — 0. 

The binding energies at typical conditions for some ele- 
ments used for Bose-condensation experiments are com- 
piled in table HI These data show that the scattering 
length can be positive and its absolute value can be of the 
magnitude or even less than ac- In these cases the bind- 
ing energy is more than two orders of magnitude higher 
than the thermal energy, i.e., the bound state is stable, 
although |(jJb| is too small compared with experimental 
values of bound states near the continuum edge. For ex- 
ample the energy of the last vibrational state below the 
continuum edge for sodium is (—13100 ± 900)neV [isl j. 
The reason for that difference is that the contact inter- 
action is a low energy approximation. One would expect 
an increase of the binding energy for an increasing poten- 
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FIG. 5: Scattering phase versus momentum for a finite-range 
two-particle interaction with one possible bound state 



have to face the paradoxical situation that bound states 
appear also for a repulsive interaction. Considering the 
BEC-BCS transition it seems therefore that one has to 
include inevitably the effect of a finite potential range in 
order to be able to describe bound states, pairing and 
BEC correctly. 
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tial range. Table U shows further that for stable bound 
states, i.e., \ujb\ ^ T, the influence of the medium on the 
binding energy is negligible. On the other hand, since a 
finite-range potential would stabilize the bound states a 
medium influence on them may be measurable. 

To summarize, it has been shown that the contact in- 
teraction produces always attractive forces, if one postu- 
lates a finite scattering length. A bound state appears 
as soon as the interaction is strong enough. Due to the 
Bose enhancement, many-body effects support the for- 
mation of bound states while thermal fluctuations and 
the motion relative to the medium hinder this forma- 
tion. The model describes the experimental experience 
that bound states and cluster formation appear before 
the Bose condensation. The calculations show that for 
bosons at finite density bound states are possible on both 
sides of a Feshbach resonance but they are quite insta- 
ble near the resonance, too. To find better agreement 
with the experiment and to make predictions for future 
experiments, the model has to be improved to describe 
an interaction with finite range and bound states with 
finite life time. Our T-matrix approximation shows that, 
in contrast to the pseudopotential ansatz, the scattering 
length is not necessarily proportional to the interaction 
strength and that whether the interaction is attractive or 
repulsive does not follow inevitably from the sign of the 
scattering length. Whether an interaction is repulsive or 
attractive can be found from the sign of the scattering 
phase, the slope of which for small momenta is related 
to the scattering length, as shown in Fig. [5] The contact 
interaction can only describe the strong and weak attrac- 
tive cases while the pseudopotential can only describe the 
weak attractive and repulsive cases. Therefore the sign of 
the scattering length alone is not decisive whether the in- 
teraction is attractive or repulsive. Otherwise one would 
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